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Abstract 

The problem of spanning trees is closely related to various interesting problems in 
the area of statistical physics, but determining the number of spanning trees in 
general networks is computationally intractable. In this paper, we perform a study 
on the enumeration of spanning trees in a specific small-world network with an 
exponential distribution of vertex degrees, which is called a Farey graph since it 
is associated with the famous Farey sequence. According to the particular network 
structure, we provide some recursive relations governing the Laplacian characteristic 
polynomials of a Farey graph and its subgraphs. Then, making use of these relations 
obtained here, we derive the exact number of spanning trees in the Farey graph, 
as well as an approximate numerical solution for the asymptotic growth constant 
characterizing the network. Finally, we compare our results with those of different 
types of networks previously investigated. 
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1 Introduction 



Counting spanning trees in networks (graphs) is a fascinating and central is- 
sue in statistical physics, because of its inherent relevance to diverse aspects 
in related fields. For example, the number of spanning trees is an important 
measure of reliability of a network [Tf2] . Again for instance, it is exactly the 
number of recurrent configurations of the Abelian sand-pile models [3f5f4"] . 
which have been studied extensively in the past two decades as a paradigm 



Email address: zhangzz@fudan.edu.cn (Zhongzhi Zhang). 



Preprint submitted to Elsevier 



23 January 2012 



of the self-organized criticality [6]T7] . On the other hand, the problem of span- 
ning trees has numerous connections with other interesting problems associ- 
ated with networks, such as dimer coverings [8], Potts model [HE], random 
walks pUfTT] . origin of fractalitity for fractal scale-free networks [T2"][T3] and 
many others [14] . 

In view of its wide range of applications, the enumeration of spanning trees 
has received considerable attentions from the scientific community. A lot of 
previous studies have focused on counting spanning trees on different me- 
dia, including regular lattices [To]ll6|ll7] . the Sierpinski gaskets [TS], and the 
Erdos-Renyi random graphs [TH], even scale-free networks [20.21] that dis- 
play the striking power-law degree distribution [22] as found for many real 
systems [23f2 4.25.26]. These researches uncovered the impacts of various ar- 
chitectures on the number of spanning trees in different networks. However, 
not all real-word networks are scale-free. It has been observed [27] that the 
degree distribution of some real-life networks (e.g., power grid) decays expo- 
nentially, although they show small-world effect [28] characterized simultane- 
ously by low average path length and high clustering coefficient. Thus far, the 
problem of spanning trees in small-world networks with an exponential degree 
distribution has not been addressed. 

In this paper, we count spanning trees in a small-world network with a con- 
nectivity distribution decaying exponentially, which is translated from the 
Farey sequence [29J and thus called a Farey graph (network). According to 
the specific structure of the Farey graph, we derive recurrence formulas for 
the Laplacian characteristic polynomials of Farey graph and its subgraphs, 
based on which we determine the exact number of spanning trees in Farey 
graph and numerical value of its asymptotic growth constant. We also com- 
pare our results with those for other networks with the same average degree 
of nodes but different degree distributions. 



2 Construction and topological properties of Farey graph 

The graph under consideration is derived from the famous Farey sequence [2H] . 
In mathematics, a Farey sequence of order n (n is a positive integer) is a 
set (denoted by F n ) of irreducible fractions between and 1 arranged in an 
increasing order, the denominators of which do not exceed n. For example, the 
first four Farey sequences are: F x = {f , ±}, F 2 = {f , \, y}, F 3 = {f , §, |, §, {}, 
-^4 = {i)I>|'2'l'I'l} - I 11 ^ ac ^' Farey sequence F n can be constructed from 
F n _i by using the Farey sum operation denoted as ©. Let | and % be two 
irreducible fractions, then one can define their "mediant" as | © ^ = 
It has been proved that F n could be obtained from F n _i by calculating the 
mediant between each pair of consecutive fractions in F n _x, keeping only those 
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mediants with denominator equal to n, and placing each mediant between the 
two values from which it was derived. The Farey sequence has an interesting 
property, that is say, any two neighboring fractions - and £ in a Farey sequence 
satisfy rq — ps = 1. 




Fig. 1. Illustration of a Farey tree for the first several levels. 

The Farey sequence can be organized in the so-called Farey tree [301131] . Be- 
ginning with the two fractions 9 and j, the first level of the tree is | = § © y, 
and the second level consists of two fractions ~ and | that are obtained by 
the Farey sum operation over all the previous fractions, i.e., ~ = ° © | and 
I = I ® r Repeating this Farey addition recursively, we can obtain the Farey 
tree, the nth level of which includes 2 n ~ l irreducible fractions, see Fig [TJ In 
addition, based on the Farey tree, one can construct a small-world graph 
called Farey graph [32|33] , in which nodes represent the irreducible fractions 
between and 1, and two nodes | and ^ are connected if they satisfy the 
relation rq — ps = 1, or equivalently if they are consecutive terms in some 
Farey sequence, see Fig[2j Notice that, the Farey tree is actually a spanning 
tree of the Farey graph. 




Fig. 2. Illustration of a Farey graph corresponding to the Farey tree depicted in Fig. 

m 

In fact, the Farey graph can be created in the following iterative way [M1I35] . 
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Fig. 3. (Color online) Alternative construction method of the Farey graph. W g +\ 
can be obtained by joining two copies of W g denoted as Wg (j) = 1,2), the initial 
nodes of which are represented by X g and Yg, respectively. In the merging process, 
Y g l and Xg are identified as the hub node Z g+ \ of W g+ \. In addition, Xg and Y g 2 
become respectively X g+ \ and Y g+ i, which are linked to each other by a new edge. 

Let W g denote the Farey graph after g iterations. For g — 0, W g is an edge 
connecting two nodes. For g > 1, W g is obtained from W g -\\ for each edge 
in Wg-\ introduced at iteration g — 1, we add a new node and attach it 
to both ends of the edge. The Farey graph is minimally 3-colorable, uniquely 
Hamiltonian, and maximally outerplanar 3(i . Particularly, the Farey graph 
exhibits some remarkable properties of real networks, it is small world with 
its average distance increasing logarithmically with its node number, and its 
clustering coefficient converges to a large constant In 2 [34~f35] . 



3 Spanning trees on Farey graph 

After introducing the network construction and its structural properties, we 
continue to investigate spanning trees in the Farey graph W g , by establishing 
and using recursive relations for the characteristic polynomials of it and its 
subgraphs at different iterations. To facilitate the description, we give the 
following definitions. Let X g and Y g denote the two nodes in graph W g that 
are created at iteration and thus called initial nodes. In addition, the node in 
W g having the highest degree was generated at iteration 1, it is thus named the 
hub node and represented by Z g . Then, the Farey graph can be alternatively 
built as follows. Given iteration g, W g+ i may be obtained by joining at the 
initial nodes two copies of W g , see Fig. [31 This new construction shows that 
the Farey graph is self-similar, implying that the network order (represented 
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by N g ) obeys N g = 2N g _ 1 - 1, which coupled with N = 2 gives N g = 2 9 + 1. 
Using the self-similar property of the Farey graph, one can enumerate the 
number of spanning trees, which is the focus of this paper. 

According to the well-known matrix-tree theorem j37j, the number of spanning 
trees in a connected graph can be expressed numerically but exactly in terms 
of the Laplacian spectra corresponding to the graph. Let Ns^(g) represent the 
number of spanning trees in the Farey network W g , whose Laplacian matrix 
is represented by L g = [L^], then one can obtain N ST (g) by computing the 
product of all non-zero eigenvalues of L g as [381139] 

1 Ng-l 

n st(9) = w Ii Hq) , (i) 

9 i=l 

where Xi(g) (i = 1, 2, . . . , N g — 1) denote the N g — 1 non-zero eigenvalues of 
matrix L g , the entry Ly of which is defined as follows: Ljj = — 1 (or 0) if nodes 
i and j are (or not) adjacent and % ^ j, equals the degree of node i ii i = j. 

Equation (CQ) shows that the issue of determining N$T(g) is reduced to finding 
the product of all nonzero Laplacian eigenvalues. Although the expression of 
Eq. ([[]) appears succinct, it requires computing the eigenvalues of a matrix 
with order N g x N g , which make heavy demands on time and computational 
resources for large networks, since the complexity for calculating eigenvalues 
is very high. Thus, for large networks, it is not acceptable to obtain N^T(g) 
through direct calculation of the Laplacian spectra, due to the limitations 
of time and computer memory. It is then of significant practical importance 
to seek for a computationally cheaper approach to overcome this problem. 
Fortunately, the self-similar architecture of the Farey graph allows for calcu- 
lating the product of its nonzero Laplacian eigenvalues to obtain an analytical 
solution to the number of spanning trees. Details will be provided below. 

We use P g (X) to denote the characteristic polynomial of the Laplacian matrix 
L g , i.e., 

P g (X) = det(L g - XI g ) , (2) 

in which I g is an N g x N g identity matrix. As mentioned above, our goal is to 
evaluate the product of all the nonzero eigenvalues of L g , namely, all nonzero 
roots of polynomial P g (X). 

In order to find the product, we denote Q g as an (N g — 1) x (N g — 1) submatrix 
of (Lg — Xl g ), which is obtained by removing from (L g — XI g ) the row and 
column corresponding to an initial node, say, node X g of network W g . In 
addition, we use H g to represent a submatrix of (L g — Xl g ) with an order 
(N g — 2) x (Ng — 2), obtained from (L g — XI g ) by removing from it the two 
rows and columns corresponding to the two initial nodes in W g , i.e., X g and Y g . 
On the other hand, let Q g (X) and R g (X) denote respectively the determinants 
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of Q g and R g . Then, the three quantities P g (X), Q g (X), and R g (\) obey the 
following relations: 



2gL - A 



^+i(A) 



-E„ 



-E s Q„ + E 5 



(3) 



and 



Q s +i(A) 



-R 9 +i(A) 



2d, 
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l q lip 

Q 3 + E 9 O 
O R„ 



2d g — A l r l r 
ij R 5 O 
L O R. f , 



(4) 



(5) 



In Eqs. (J3KSJ) , E 9 is an (N g — 1) x (N g — 1) square matrix with only one entry 
equal to 1 while all other entries being 0, viz., 



E n 



1 ••• (A 
••• 



\0 •■• 0/ 



(6) 



2d 9 expresses the degree of the hub node Z g+ i in the (g + l)th generation 
W g+ i, which is in fact the double of the degree, d g , of an initial node X g or Y g 
in graph W g ; l q (l r ) is a vector of order N g — 1 (N g — 2) with d 9 entries being 
— 1 all other N g — d g — 1 (iV 9 — d g — 2) entries equaling 0, in which each entry 
— 1 describes an edge connecting the hub node Z g+ i and a node belonging 
to either W g or W g , both of which are amalgamated into finally, the 

superscript T of a vector represents transpose. 

For the convenience of following description, we index respectively the rows 
and columns of P g+ i(A) (or its variants obtained from it by elementary matrix 
operations) by Vi and C« with i = 1, 2, ... , 2N g — 1. Applying the elementary 
matrix operation, i.e., replacing column Cj(N g + 1 < j < 2N g — 1) by Cj — 
Cj-N g +i, we have 
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Pg+lW 



2d g - A 

c 



i 



o 



Q9 + — (Q 9 + 2E fi 

"Eg Qg + 2E fl 



(7) 



We continue to replace the row V%{2 < i < N g ) in Eq. (J7|) by Vi + Vi + N g -i, 
yielding 



2d g - A l q 



2ZJ 



o 



2c?g — A 



2/ T 



Eg Qg + 2Eg 

Qg + 2E 9 | 



Thus, we have expressed the P g+ i(X) as a product of two determinants, which 

we denote by P fl +\(A) and Pg+i(\), respectively. These two determinants can 
be easily evaluated as 



^i(A) 



2dg — 2A lq 


+ 


A 


211 §9 




2/J Q 3 



2P g (A) + Ag,(A) 



(9) 



and 



Pg+lW = \Qg\ + 2 |R S | = Qg(\) + 2R g (X). 



(10) 



Inserting Eqs. (Q and (ITOj) into Eq. ([8]), we have 



P g+l (X) = [2P,(A) + Ag 9 (A)][g 9 (A) + 2P 9 (A)] 
= 2P 9 (A)g 9 (A)+4P 5 (A)P,(A) 
+A[(g,(A)) 2 + 2g,(A)P 9 (A)]. 



Similar to the computation processes for P g +i(X) and P 9 +i(A), g g +i(A) and 
R g+ i(X) can be calculated as shown in Eqs. (jT2"|) and (TT5|) . 
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d„ — A 



In 







Qg + E g O 

O R„ 



+ 



A 

q Q 9 + E, O 



c 



O R„ 



d n - A 



Z r 
Q 9 + E 9 O 
O R„ 



R 9 (X)[P 9 (X) + Q g (X)} + Q g (A)[Q s (A) + R g (X)) 

+XRg(X)[Qg(X) + Rg(X)], 



:i2) 





dg - 


A l r 







dg ~ 


A 






A 





Rg+iW - 


17 


Kg 


O 






Kg 


O 




;T r> 

r ^g 









o 


Rg 




17 


o 


Rg 




17 


Rg 



= Q g (X)R g (X) + Q g (X)R g (X) + X[R g (X)] 2 . (13) 

Having derived the recursive relations for P g (X), Q g (X), and R g (X), shown 
in Eqs. (II lti!3j) . we proceed to compute the product of the nonzero roots of 
polynomial P g (X). Since P g (X) has one and only one root equal to zero, say 
Ao(<?) = 0, to find this product, we define a new polynomial P g (X) as 

^(A) = \Pg{X) . (14) 

Then, it is evident that 

Ng—l Ng-1 _ 

II Hg)= II Hg), ( 15 ) 

i=l i=l 

in which Xi(g), X~2(g), ■ ■ ■ , X~N g -i(g) represent the N g — 1 roots of polynomial 
P g (X). Thus, the determination of the product of nonzero eigenvalues of Lapla- 
cian matrix L 5 is equivalent to calculating the product on the right-hand side 
(rhs) of Eq. (ifej). 

To find the product Il^i 1 Aj(g), we express polynomial P g (X) in the following 
form, i.e., P g (X) = Ylflo 1 Pgij)^ ■> XTl which p g (j) is the coefficient of term A J 
with degree j. Since it is obvious that p g {N g — 1) = —1, we then have 

Ng-l Ng-1 

e p 9 m j = - n t A - • (is) 

3=0 i=i 
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According to Vieta's formulas, the following relation holds: 

Ng—l _ 

I] Mg) = -PM- (17) 



t=i 



Thus, all we need is to determine the constant term p g (0) of polynomial P g (X). 
From Eqs. (I11H13P it is not difficult to derive the following recursion equations: 

P g+1 (X) = 2P g (X)Q g (X)+AP g (X)R g (X) 

+ [Q g (X)] 2 + 2Q g (X)R g (X), (18) 

Qg+iW = [Q 9 W? + 2Q s (A)i?,(A) + X[P g (X)R g (X) 

+Q g (X)R g (X) + (R g (X)n (19) 

and 

R 9+1 (X) = 2Q g (X)R g (X) + X[R g (X)] 2 . (20) 

On the basis of above relations, we can find the value for p g (0). To this end, 
we give some additional variables. Let q g (0) and r g (0) be the constant terms of 
Q g (X) and R g (X), respectively. According to Eqs. fll8H20p . the three quantities 
p g (0), q g (0) and r g (0) obey the recursive relations: 

p g+1 (0) = 2p g (0)q g (0) + # 9 (0)r 9 (0) + [q g (0)} 2 + 2g,(0)r 9 (0) , (21) 

q g+1 (0) = [q g (0)] 2 + 2q g (0)r g (0), (22) 

and 

r g+1 (0) = 2g 9 (0)r,(0). (23) 

Plugging Eq. fl23J into Eq. (|22j) to obtain 

g s+1 (0) = [g 9 (0)] 2 + r 9+1 (0), (24) 
which can be rephrased as 

r s+ i(0) = g 5+ i(0)-[g 9 (0)] 2 . (25) 

Replacing r g (0) in Eq. (1231) by the expression given on the rhs of Eq. (1251) 
leads to 

q g+1 (0) = 3[q g (0)] 2 - 2g,(0)[g 3 _ 1 (0)] 2 . (26) 



9 



Thus, we obtain the recursive relation governing q g+ i(Q), q g (0), and g 9 _i(0), 
as shown explicitly in Eq. (|26|) . 

Solving Eq. (j2BJ) one can arrive at the formula for q g (0). For this purpose, we 
introduce an intermediary quantity k g , defined as 

k g = q a (0)/{q a - 1 (0)}\ (27) 
making use of which Eq. ( 1261) can be rewritten as 

K+i = 3 - £ . (28) 

Considering the initial condition k\ = gi(0)/[go(0)] 2 = 3, Eq. ff28l can be 
solved to yield 

k = 2(-l)'-(-2)-' (29) 

9 (-1)9 - {-2)~9 ■ 1 J) 

With the obtained exact result for k g , we can reword Eq. ( 1271) as 

\nq g (0) = 21ng 9 _!(0) + lnk g . (30) 

Using the initial condition lngx(0) = In 3 and the expression for k g provided 
by Eq. ( 1291) . Eq. ( 1301) can be solved inductively to obtain 

\nq a (0) = 2 9 " 1 E2^ In . (31) 

i=o z l 

Thus, we have 



3-1 /oi+2 _ 1 \ 

5.(0) =n 



i=0 



2 i+l 



1 

9 



(2^ +1 - l) fl (2 J - l) 



29-' 



(32) 



i=2 



After deriving g 9 (0), we now are in position to calculate p g (0). Notice that 
Eq. (12T|) can be decomposed into a product of two terms as 

p g+1 (0) = [2 Pg (0) + q g (0)][q g (0) + 2r g (0)} . (33) 

In addition, Eq. f[2"2"j) can be rewritten as 

q g+1 (0) = q g (0)[q g (0) + 2r g (0)]. (34) 
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Then, we have 

fc+i(0) q g (0) 1 J 

which together with Pi(0)/gi(0) = —3 leads to 

H = -(2' + D. (36) 
Thus, we can obtain the explicit formula for p g (0): 

p g (0) = - (2* + 1) (2^ - l) f[ (? - I) 29 "' . (37) 

i=2 



Hence, the number of spanning trees in the Farey graph W g is: 

Nst(9) = = - l) f[ (2 l " if" • (38) 

Equation (|38|) is our main result, which is exact and holds for any legal g. 
Note that Eq. (|38|) can also be derived using another approach provided in the 
appendix. We thank an anonymous referee reminding us of this nice method. 

Using the above-obtained result given by Eq. (1381 . one can determine the 
asymptotic growth constant of the spanning trees — an important quantity 
characterizing network structure — for the Farey graph, which is defined as 
the limiting value |4"0f4"T] 

E Wa = lim = lim (39) 

that converges to a constant value 0.9458, a finite number a little smaller than 
1. 

The obtained entropy of spanning trees for W g can be compared to those 
previously found for other media with the same average node degree as the 
Farey network. In the pseudofractal fractal web, the entropy is 0.8959 [20], a 
value less than 0.9458. For the square lattice and the two-dimensional Sier- 
pinski gasket, their entropy of spanning trees are 1.16624 [15] and 1.0486 [18], 
respectively, both of which are greater than 0.9458. Therefore, the number 
of spanning trees in the Farey graph is larger than that of the pseudofractal 
fractal web, but is smaller than that corresponding to the square lattice or the 
two-dimensional Sierpinski gasket. The distinctness lies with the architecture 
of these networks. Although they have identical average node degree, they 
show quite different degree distributions. Thus they exhibit disparate distri- 
bution of Laplacian spectra that have been shown to display similar distribu- 
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tion as node degrees [4~2j43f44j . which fundamentally determine the number 
of spanning trees. 



4 Conclusions 

In this paper, we have investigated the problem of spanning trees in a Farey 
graph with the small-world effect and an exponential degree distribution. On 
the basis of its structure self-similarity and a decimation procedure, we de- 
rived some recursion relations of the Laplacian characteristic polynomials for 
the Farey graph and its subgraphs at different iterations. We then applied 
these useful relations to enumerate spanning trees in the Farey graph and ob- 
tained the exact number of spanning trees, as well as the numerical result of 
asymptotic growth constant. An advantage of our technique lies in the avoid- 
ance of laborious computation of Laplacian spectra that is needed for a generic 
method for determining spanning trees in general networks. We also compared 
our results with those previously obtained for other networks. Finally, it should 
be mentioned that there are many interesting questions about the Farey graph 
for future research, e.g., determining the number of its sub-trees when only 
nodes with denominator less than a given value n are kept. 
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A An alternative method for determining the number of spanning 
trees in the Farey graph 

Here we introduce simply the idea of another method for enumerating span- 
ning spanning trees in graph W g . Suppose that one considers a spring hamil- 
tonian on W g , where all spring constants are equal, and equal to K. Then the 
spring Hamiltonian is 



in which Xj is the scalar displacement at node % of the graph, and sum extends 
over all edges i ~ j of the graph. 




(A.l) 
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If we calculate the partition function of this graph, by integrating over all 
Xi, except one node is kept at zero displacement, the partition function is 
easily seen to involve the determinant of the corresponding Laplacian matrix. 
Suppose, we integrate over all nodes except the two initial nodes. Then, the 
restricted partition function can be written in the form 

W g (x,y) = A g e W (-K g {x-y) 2 ) , (A.2) 

where x and y are the displacements at the two initial nodes, and K n may be 
called the effective spring constant between them. 

The self-similar structure of the Farey graph implies that we can express A g+ \ 
and K g+ i in terms of A g and K g . In fact, 

W g+1 (x,y) = J dzW g (x,z)W g (z,y), (A.3) 

where z is the displacement attached to the hub node. From Eq. (|A.3|) . it is 
easily seen that K g+ \ = \K g + 1. This relation is easily solved explicitly to 
obtain K g , and we can put this solution into the recursion equation for A g to 
get Eq. f[3"51) in the main text. 
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